Ab initio description of atomic quantum gases interacting with a single carbon 

nanotube 
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We study inelastic processes in the hybrid quantum system constituted by a carbon nanotube 
(CNT) in contact with an ultracold quantum gas, such as a cloud of thermal atoms or a Bose- 
Einstein condensate (BEC). We present a parameter- free ab-initio approach for the loss rate based 
on the underlying scattering process, considering the two-dimensional character of the system as 
well as the exact Casimir-Polder potential. The predicted loss rates are in perfect agreement with 
recent experimental results, obtained both for a thermal cloud of rubidium atoms and for a BEC. For 
the trap loss of a thermal cloud, we find that retardation effects become important and contribute 
significantly, which emphasises the crucial role of the exact interaction potential. 

PACS numbers: 78.67.Ch, 67.85.-d, 34.50.-s, 82.20.Ej 
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During the last years advanced techniques of cooling 
and controlling single atoms or ions, as well as clouds 
or condensates, have successfully been combined with 
sophisticated solid-state techniques in the production 
of mesoscopic structures on the nanometer scale. The 
emerging variety of hybrid quantum systems constitutes 
an interesting new field of physics [l|-|j]. Nanotubes, 
nanowires and nanorods have proven effective for con- 
structing such hybrid systems [5j-l7|. These can be com- 
bined to carpets of dense standing nanotubes forming a 
structured surface 0, 0] or to an array of nanotubes serv- 
ing as a diffraction grating 13, 11 1. Particular attention 
is currently being given to the fundamental system of a 
single nanotube interacting with a cold gas of atoms 
[lH | , which lays the foundation for understanding systems 
of higher complexity. This is by far not a trivial problem 
and raises the question of the applicability of established 
scattering theory and non-dynamical approaches (ljj ). 

In a recent experiment [l2l . Il7j Schneeweiss et al. mea- 
sured the losses of atoms absorbed by a multiwall carbon 
nanotube (length L = 10.25/im and diameter ranging 
from 275nm at the bottom to 40nm at the tip) grown 
on top of a nanochip. The CNT was immersed into a 
BEC and into a cloud of thermal atoms and the trap loss 
7 was measured as a function of the distance d between 
the center of the trap and the surface of the nanochip 
(see Fig. [T]). 

In this Letter we present a parameter-free ab-initio 
description for the absorption of ultracold atoms by 
cylindrical geometries, based on the underlying scatter- 
ing process. We show that an accurate description is 
possible, if the exact interaction potential is considered 
and the scattering process is treated properly. However, 
this remains non-trivial in the present case, due to the 
two-dimensional character of the scattering process 0- 
[20l | and due to the intricacy of the exact interaction 
potential 14|, 1 2 1 1 . In contrast to fitting an arbitrarily 
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FIG. 1. (Color online) Schematic view of a CNT immersed 
into a cloud of trapped atoms. The nanotube is grown on top 
of a nanochip and the center of the cloud is held at a distance 
d to this surface. 



chosen model potential to the experimantal data [121 Il7| , 
this parameter-free approach offers a foundation for 
extensions to more complex hybrid systems and might 
assist the design of CNT based nanodevices. 



Descriptions of an atom interacting with a cylindrical 
geometry go back to Zel'dovich [22| who analyzed the 
interaction between an atom and a perfectly conducting 
cylinder in 1935. Further work on this subject extended 
the potential to the the more realistic case of a cylin- 
der with finite conductivity. A closed form of the full 
Casimir-Polder potential between an atom and a dielec- 
tric cylinder of radius R was first given by Nabutovski 
et al. in 1979 [23| and later by Marvin and Toigo in 
1982 [24|]. Though this problem has been revisited over 
the last years by several groups 25 27], the result ob- 
tained by Nabutovski et al. [23j is still the most general, 
as it covers the full Casimir-Polder potential of an atom 
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interacting with a cylindrical geometry with arbitrary 
dielectric properties, including the perfectly conducting 
case. 

The Casimir-Polder potential for a dielectric cylinder 
exhibits a nontrivial transition from the well-studied 
atom-wall potential at small distances (r — R) from 
the atom to the surface of the cylinder, which is 
— C-3/(r — R) 3 , to the fully retarded asymptotic —Cg/r 6 
behavior far away from the cylinder. The coefficients C3 
and Cq, are well defined quantities depending on both 
the dielectric constant of the tube and the properties of 
the atom, with Cr additionally depending on the radius 
R of the cylinder [23j . The full Casimir-Polder potential 
is approximated by the van der Waals potential only 
at very small distances and might resemble the van 
der Waals long-range asymptote — Cs/r 5 only in a very 
narrow transition region, if at all. An accurate method 
for the numerical treatment of the full Casimir-Polder 
potential and, in particular for the nontrivial transition 
between the van der Waals and highly retarded limit, 
has only been developed very recently 13L 21|. 



The interaction potential between an atom and an ac- 
tual CNT of finite length L obviously differs from the 
atom-cylinder potential fjjj in a nontrivial fashion and 
remains i n g eneral unknown; pairwise-summation ap- 
proaches |28| can not reproduce the correct Casimir- 
Polder forces [2^. Variations of the interaction energy 
along the tube axis are, however, small as long as the ra- 
dius R(z) of the CNT varies smoothly. In this case, the 
longitudinal free motion of an atom along the CNT may 
be separated from the two-dimensional dynamics perpen- 
dicular to the tube, which is governed by the interaction 
potential 23| of the atom with a cylinder of radius R{z). 
The contribution of a single perpendicular plane to the 
total loss rate 7 can be given in a differential form 



d 7 (z) = n(0, 0, z)K™& z 



(1) 



with the density of particles n around the nanotubc, lo- 
cated at (0,0) and the loss rate constant Kf^ for inelas- 
tic and reactive scattering in this two-dimensional (2D) 
subsystem. These inelastic reactions (sticking, adsorp- 
tion, ...) happen at short distances to the surface and are 
described in the Langevin model where all atoms reaching 
the surface contribute to the rate constant K?® [30|, [3l[ . 
The total trap loss 7 of the full system is obtained by 
integration over all 2D planes, 
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n(0,0,z)K™dz. 



(2) 



Deviations from this description might be expected at 
the tip of the nanotube where on the one hand, the 
dynamics along the tube axis is influenced by the atom- 
cylinder potential and on the other hand, the potential 
differs from the infinite cylinder potential. However, 
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FIG. 2. (Color online) Trap loss 7 for a BEC of rubidium 
atoms overlapping with a carbon nanotube. The calculated 
trap-loss for a dielectric constant e = 5 (solid line) reproduces 
the data given in [l2| | (full circles and errorbars). The dot- 
dashed (dashed) line shows results for a dielectric constant of 
e = 2.5 (e = 100) which leads to rather small variations in 
the predicted loss rates, compared to the experimental error 
bars. 



these deviations will only give small corrections to the 
total trap loss. 

Within this description, we are able to calculate the 
trap loss of a Bose-Einstcin condensate (BEC) and of 
a thermal cloud of atoms interacting with a CNT, as 
measured for both cases by Schneeweiss et al. [l2[. All 
calculations presented below are based on the exact 
experimental parameters 
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of the CNT radius R(z) is assumed to be linear. The 
center of the cloud is shifted away from the tube in 
y-direction and lies on a ray forming an angle of 13° 
with the z-axis [HI ■ Furthermore, the experimental data 
is, for both cases (BEC and thermal cloud), assigned 
to trap-surface distances measured separately with a 
BEC (HHj. 

Figure [2] shows the trap loss 7 for a BEC of rubidium 
atoms as a function of the trap-surface distance d. The 
BEC is described by a macroscopic wave function in the 
ground state of the trap and the presence of an absorb- 
ing impurity in the cloud leads to a local depletion of 
its density. As a result, the condensate readjusts with 
the characteristic speed of sound v s (here 0.8mm/s) j3(| 
which causes a flux toward the nanotube [StJ • Due to the 
symmetry in the region close to the nanotube this flux is 
symmetric around the tube axis. Only s-waves (m = 0) 
reflect this symmetry and contribute to the absorption 
cross section <7 a b s - An additional factor 27r accounting 
for the isotropy of the incoming flux has to be considered 
and the loss rate constant is 



t^2D o (m=0) / n 



(3) 



The density of the BEC in a harmonic trap is given by a 
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FIG. 3. (Color online) Trap loss 7 for a thermal cloud of ru- 
bidium atoms overlapping with a carbon nanotube. The black 
solid line shows the predicted trap loss for a dielectric constant 
of e = 5 which is in good agreement with the experimental re- 
sults [T2| (full circles and errorbars). The blue triangles show 
a three-dimensional classical simulation of the system which 
reproduces the predicted classical trap loss (blue dashed line) 
and shows rather small deviations from the quantum mechan- 
ical results. 
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FIG. 4. (Color online) Absorption cross section for rubidium 
atoms scattered at a cylinder of radius 50nm and a dielec- 
tric constant of £ = 5. The black curve shows the full ab- 
sorption cross section <7 a b s including all contributing partial 
waves. The blue and green curves show the individual contri- 
bution of each partial wave from m = (s-wave) up to m = 4 
(g-wave). The orange region in the background shows the 
velocity distribution perpendicular to the tube (in arbitrary 
units) for a cloud of lOOnK. Notice that due to the reduced 
dimensionality of the system <r a b s is a length. 



Thomas-Fermi distribution [36[, with the characteristic 
Thomas-Fermi radii, rJ F = 16/im, F = rJ F = 3.2/im. 
As the number of atoms in the condensate decreases 
during the absorption process, the Thomas-Fermi radii 
and the speed of sound do not stay constant over time, 
which leads to an algebraic instead of an exponential 
decay [37[. These corrections and possible fluctuations 
due to collective oscillations [12], 15 1 lie within the 
experimental error bars. The influence of the nanochip 
on the trap loss can be neglected for distances d larger 
than 5^m 



12j . The cross section a 



(m=0) 



is calculated 

using the full Casimir-Poldcr potential [231 ] and incoming 
boundary conditions for the description of inelastic 
collisions at the surface of the nanotube [3, HH . The 
experimental data [l2[ shown in Fig. [2] (full circles) is 
perfectly reproduced by our calculations for a dielectric 
constant e — 5 (solid line). Variations of the dielectric 
constant in the range of e = 2.5 (dot-dashed line) up to 
e = 100 (dashed line) only lead to small changes in the 
calculated loss rates. 

The trap loss for a thermal cloud of rubidium atoms is 
shown in Fig. [3] The velocity of the atoms in a thermal 
cloud is given by a Maxwell-Boltzmann distribution with 
T = lOOnK and the loss rate constant is given by 



K 2D = 



V±cr abs (t>jJ) 



(4) 



where v± is the velocity in the plane perpendicular to the 
nanotube. The density of the atoms in the thermal cloud 
is given by a gaussian distribution with a x = 27 fim and 
o y = cr z = 5.5/im. 

The absorption cross section er a b s in two dimensions 
for rubidium atoms scattered by an infinite cylinder of 



radius 50nm and dielectric constant e = 5 is illustrated 
in Fig.|4]and has been calculated using the full interaction 
potential [23| and incoming boundary conditions (l4. 38. 
|39|. For a cloud with T = lOOnK, angular momenta up 
to the (/-wave (m = 4) contribute significantly to the 
absorption cross section er a b s at a CNT with a radius of 
50nm, which is in contrast to simpler models [l2j . 

Our prediction shown in Fig. [3] (solid line), which is 
based on the full Casimir-Poldcr potential and a dielec- 
tric constant e = 5, is in good agreement with the ex- 
perimental results (full circles). Variations of e lead to 
rather small variations of our results compared to the 
experimental error bars (see Fig. [5]). 

For the case of a thermal cloud we also performed a 
time dependent three-dimensional classical simulation 
of a cloud of 1000 atoms at constant temperature in a 
harmonic trap, overlapping with an absorbing nanotube. 
The interaction potential between the atoms and the 
tube is approximated by the long-range asymptotics of 
the Casimir-Poldcr potential. The effect of the finite 
length of the tube has been taken into account by 
a smooth decrease of the potential at the endings of 
the tube. The loss rate observed in this simulation 
(triangles in Fig. [3]) is reproduced by the result predicted 
via Eqs. ^ and (@| using the corresponding classical 
absorption cross section for the asymptotic potential 
(dashed line in Fig. [3]). The comparison of these results 
with the full quantum mechanical calculation and the 
experimental data shows a good agreement. Thus, in 
contrast to the extremely nonclassical BEC [see Eq. (|3J)], 
the behaviour of the thermal cloud at T = lOOnK is 
essentially determined by classical dynamics. 
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FIG. 5. (Color online) Trap loss 7 for a thermal cloud of ru- 
bidium atoms. The loss rates has been calculated via Eq. |[2| 
based on the full Casimir-Polder potential (lower curves) and 
on the van der Waals potential only, where retardation ef- 
fects are neglected (upper curves). The shaded area shows 
loss rates for dielectric constant of the nanotube in a regime 
from e = 2.5 (dot-dashed line) up to e — 100 (dashed line). 
The solid line shows the trap loss for a dielectric constant of 
e — 5. The experimental data [12] (full circles) clearly devi- 
ates from the prediction calculated with the van der Waals 
potential only. 



A fundamental question that needs to be addressed 
in future experiments is the influence of retardation on 
scattering processes for such systems. It has recently 
been shown that for realistic atom-nanotube systems, 
retardation effects are important in the case of a 
perfectly conducting tube [13j . This result holds for 
the present case of a carbon-nanotube immersed into a 
cloud of thermal atoms; the nonretarded van der Waals 
potential fails to describe the trap loss (upper curves 
in Fig. [5J, within a wide range of dielectric constants 
from e = 2.5 (dot-dashed line) up to e = 100 (dashed 
line). The results obtained with the full Casimir-Polder 
potential (lower curves in Fig. [5]) match the experimental 
data within the same range of dielectric constants. For 
a BEC the differences of the loss rates obtained from 
the van der Waals potential and the full Casimir-Polder 
potential are rather small and not resolvable in this 
range of dielectric properties. 

In this Letter, an accurate parameter-free ab-initio 
description of the trap loss in a hybrid system consisting 
of a single nanotube immersed into an ultracold atomic 
quantum gas, such as a thermal cloud or a BEC is 
presented. A quantum mechanical calculation of the ab- 
sorption cross section based on the exact Casimir-Polder 
potential leads to perfect agreement with recent exper- 
121 ] . both for a thermal cloud and for a BEC. 



in the Casimir-Polder potential need to be considered. 
An accurate description of this hybrid system is achieved 
with a proper use of scattering theory together with the 
exact Casimir-Polder potential. The present approach 
gives insight into the underlying processes and therefore 
promotes a deeper understanding of hybrid systems 
which is essential for the design of future nanodevices. 



The authors thank J. Fortagh, A. Giinthcr, P. 
Schnceweiss, T.E. Judd for providing experimental pa- 
rameters and for helpful discussions. Useful input by 
H. Friedrich is gratefully acknowledged. This work was 
supported by the Deutsche Forschungsgemeinschaft, Az.: 
FR591/13-2. 



imcnts 

Furthermore, it has been shown that the van der Waals 
potential fails to reproduce the loss rates of a thermal 
cloud of atoms; retardation effects that are accounted for 
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